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ABSTRACT 



The early-type spiral NGC 4698 is known to host a nuclear disc of gas and stars which is 
rotating perpendicularly with respect to the galaxy main disc. In addition, the bulge and main 
disc are characterised by a remarkable geometrical decoupling. Indeed they appear elongated 
orthogonally to each other. In this work the complex structure of the galaxy is investigated 
by a detailed photometric decomposition of optical and near-infrared images. The intrinsic 
shape of the bulge was constrained from its apparent ellipticity, its twist angle with respect 
to the major axis of the main disc, and the inclination of the main disc. The bulge is actually 
elongated perpendicular to the main disc and it is equally likely to be triaxial or axisymmet- 
ric. The central surface brightness, scalelength, inclination, and position angle of the nuclear 
disc were derived by assuming it is infinitesimally thin and exponential. Its size, orientation, 
and location do not depend on the observed passband. These findings support a scenario in 
which the nuclear disc is the end result of the acquisition of external gas by the pre-existing 
triaxial bulge on the principal plane perpendicular to its shortest axis and perpendicular to the 
galaxy main disc. The subsequent star formation either occurred homogeneously all over the 
extension of the nuclear disc or through an inside-out process that ended more than 5 Gyr ago. 

Key words: galaxies: bulges - galaxies: formation - galaxies: individual: NGC 4698 - galax- 
ies: photometry - galaxies: spiral - galaxies: structure. 



1 INTRODUCTION described by dynamical models of oblate ellipsoid s which are flat- 

tened by rotation with little or no anisotropy (e.g., |pTgnatelli et al.l 

Several mechanisms nave been proposed for bulge assembly. u^.i rr r. — : — T~r\U™A • r i ■ i . 

r r a J 2001; Cappellari et al. 2006), the twisting of the bulge isophotes 

Bulges formed either before the disc by hierarchical merging, or ^ blad " 1956; Zarltsky & Lo 1986), the presence of non- circular 

at the same time in a monolithic collapse, or after the disc as a , |„ . ,lr. nn J K- i - — 7; — — , l™r,A 

, , _ , , . ., , gas motions (e.g., Coccato et al. 2004; Falcon-Barroso et al. 2006), 

results of secular evolution. Furthermore, there is evidence that ... ... . , . .. . r , , , 

and the misalignment between the major axes of the bulge and 
bulges can experience accretion events via infall of external ma- disc ( Be rt ol a , Vietri7& Zelh^[l99lt lillendez-Abreu et afco! 
terial or satellite galaxies. None of these scenarios alone is able , „ ... „„. , , . , , . , , . , 

hereafter MA+08) observed in several galaxies cannot be explained 
to reproduce the observed properties of all the bulges along the . r , , ... , t , . . _. . 

r r r a ° it the bulge and disc are both axisymmetric. These features are 

Hubble sequence since bulges are a class of diverse and heteroge- . . . , . r , , . ... „ - . 

- interpreted as the signature of bulge tnaxiahty. Perfect axisym- 

neous obiects. The large bulges of lenticulars and early-type spirals . , , ..... . . . . , 

metry is ruled out when the intrinsic shape of bulges is deter- 
are similar to low-luminosity elliptical galaxies, whereas the small , , . ., .. . , . , , , . , , ... .. 

J r c mined by statistical analyses based on their observed elliptic- 

bulges of late-type spirals are reminiscent of disc components (see jFathi & Peletied 2003, MA+08). A large fraction of bulges 

Kormendy & Kenmcutt 2004 and references therein). , , , ... . , . , t . , 

: — — : ~ . — : — ! are characterised by an elliptical equatorial cross-section, and 

The study of the intrinsic shape of bulges is a key constraint „ .. „ ' , , , , ... 

J raj most of them are flattened along the polar axis they share with 

on their formation. Although the kinematics of many bulges is well ,. ,. . " . . . . , . , 

° the surrounding disc. As a consequence, triaxial bulges elon- 

gated per pendicularly with respect to the disc are expected to be 



gated perpendicularly with respect to the disc are expected to be 
very rare (rMendez-A breu et all |20 1 Ol hereafter MA+10). To date 
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NGC 4698 dBertola et aljl999h. NGC 46 72 dSarzi et alfcooob . and 
UGC 10043 dMatthews & de Griisll2004l) are the only spiral galax- 
ies known to host a prominent bulge sticking out from the plane of 
the disc. 

We decided to revisit the case of NGC 4698 since such a 
rare galaxy represents an excellent test case to derive the intrinsic 
shape of the bulge by applying the method of M A+10. NGC 4698 
is an Sab(s) spiral dde Vaucouleurs et al.lll99ll hereafter RC3) in 
the Virgo Cluster. Its total B-band magnitude is Bt = 11.46 
mag (RC3), which corresp onds to Mb = —19 .69 mag assum- 
ing a distance of 17 Mpc (Freedma n et alj[l994l) . The inner re- 
gion of the large elliptical-like bulge is elongated perpendicular 
to the major axis of the main disc, i.e. it is a polar bulge; this 
is also true for the outer parts of the bulge if a parame tric pho- 
tometric decomposition is adopted (Bertola et al. 1999). At the 
same time, t he stellar dBertola et al.ll 19991 : ICorsini et al.lll999T) and 
ionised-gas dBertola & Corsinil2000h components are characterised 
by an inner velocity gradient and a central zero-velocity plateau 
along the minor and major axes of the disc, respectively. This 
kinematically-decoupled core is rotating perpendicularly with re- 
spect to the galaxy ma in disc. It correspo nds to a nuclear stel- 
lar disc (hereafter NSD; |Pizzella et all 2002). The geometrical and 
kinematical orthogonal decoupling of NGC 4698 can hardly be ex- 
plained without invokin g the acquisition of exter nal material from 
the galaxy outskirts (see lBertola & Co rsini 2000). 

In this Letter, we improve the previous results by quantita- 
tively constraining the intrinsic shape of the bulge of NGC 4698 
in Section[2] The analysis of a near-infrared (NIR) image from the 
United Kingdom Infrared Telescope (UKIRT) allows us to map the 
mass distribution and minimise the contamination from dust. We 
get a more reli able estimate of the galaxy's structural parameters 
with respect to iBertola et al.l d 1999b since they folded the eastern 
side of their optical image around the galaxy major axis to deal 
with the strongest dust lanes of the main disc. In addition, they 
assumed the galaxy surface brightness to be the sum of a de Vau- 
couleurs bulge and an exponential disc with orthogonal major axes. 
Here, no a priori choice is made about the orientation of th e bulge. 
The structure of the NSD detected bv lPizzella et al.ld2002l) is anal- 
ysed in Section[3]in greater detail. The photometric decomposition 
of multi-band optical images of the galaxy nucleus obtained with 
the Hubble Space Telescope (HST) allows us to investigate the for- 
mation process of the NSD by studying its stellar populations. Our 
findings are discussed in Section|4] 



2 THE BULGE 

2.1 United Kingdom Infrared Telescope imaging 

A K— band image of NGC 4698 is available in the first data re- 
lease of the La rge Area Survey of UK IRT Infrared Deep Sky Sur- 
vey (UKIDSS. .Lawrence et alj|2007f) . It was obtained on 23 Jan- 
uary 2007 with a total exposure time of 40 s. The 3.8-m UKIRT 
telescope is operated in Mauna Kea Observatory (Hawaii, USA). 
It mounted the Wide Field Camera (WFCAM) with four Rock- 
well Hawaii-II devices. Each of them consists of 2048 x 2048 
HgCdTe detectors of 18 x 18 fim 2 and covers a field of view of 
13.7 x 13.7 arcmin 2 with an image scale of 0.4 arcsec pixel -1 . 
The gain and read-out noise were 4.5 e~ count -1 and 25 e~ 
rms, respectively. The image was reduced, cleaned for cosmic 
rays, sky subtracted, and flux calibrated with the WFCAM pipeline 
(Ham blvetal]|2008l) . A two-dimensional fit with a circular Gaus- 
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Figure 1. Two-dimensional photometric decomposition of NGC 4698. Top 
row: UKIDSS i<"-band image (left-hand panel), best-fit image (central 
panel), and residual (i.e., observed-model) image (right-hand panel). North 
is up and East is left. Bottom row: Ellipse-averaged radial profiles of the 
surface brightness (left-hand panel), ellipticity (central panel), and position 
angle (right-hand panel) measured in the UKIDSS (dots) and model im- 
age (green continuous line). The dashed blue and dotted red lines represent 
the intrinsic surface-brightness radial profiles of the bulge and main disc, 
respectively. 



sian to the field stars in the resulting image yielded a FWHM = 
0.76 arcsec. 

The photometric decomposition of the WFCAM image of 
NGC 4698 was performed using the Galaxy Surface Photom- 
etry Two-Dimensional Decomposition (GASP2D) algorithm by 
MA+08, which yields the structural parameters for a Sersic bulge 
and an exponential main disc. 

The fitting algorithm relies on a x minimisation of the in- 
tensities in counts, for which we must adopt initial trial parameters 
that are as close as possible to their final values. The trial values 
of the effective surface brightness, effective radius, shape param- 
eter, position angle, and axial ratio of the bulge (I c , r c , n, PAb, 
and gb) and of the central surface brightness, scalelength, position 
angle, and axial ratio of the main disc (Io, h, PA<j, and q^) were 
estimated from the analysis of the ellipse- averaged radial profiles 
of surface brightness /j,k, ellipticity e, and position angle PA. The 
latter were measured using the IRAF task ELLIPSE and analysed 
by following the prescriptions of MA+08. Starting from these ini- 
tial trial parameters the different photometric models of the surface 
brightness were fitted iteratively to the galaxy image. Each image 
pixel was weighted according to the variance of its total observed 
photon counts due to the contribution of both galaxy and sky, and 
determined assuming photon noise limitation and taking the detec- 
tor read-out noise into account. Seeing effects were also taken into 
account by convolving the model image with a circular Gaussian 
point spread function (PSF) with a FWHM matching the observed 
one. The convolution was performed as a product in Fourier do- 
main before the least-squares minimisation. Only the image pixels 
with an intensity larger than 0.5 times the sky standard deviation 
were included in the fit. Foreground stars were masked and ex- 
cluded from the fit. The result of the photometric decomposition is 
shown in Fig.Q] The residual image shows the tightly-wound spiral 
arms aligned with the disc major axis. Nevertheless, the fit was sat- 
isfactory. The model surface brightness is consistent with the data 
within the error bars. The differences of ellipticity and position an- 
gle between the observed and model isophotes are smaller than 0. 1 



©2012 RAS, MNRAS 000, iTfJ-lJg] 



NGC 4698 3 



and 3°, respectively. The PA is scattered in the radial range between 
about 5 and 20 arcsec where the isophotes are nearly round. 

The errors on the fitted parameters were estimated through 
Monte Carlo simulations on images of artificial galaxies. A set of 
disc galaxies with a Sersic bulge and an exponential main disc with 
7 < K T < 8 mag to bracket the magnitud e of NGC 4698 (K T = 
7.54 mag. lKassin. de Jong. & Poggdl200q) was generated. The ar- 
tificial galaxies were assumed to be observed at the distance of the 
Virgo Cluster taking into account resolution effects. The parameters 
of the artificial galaxies were rando mly chosen in the ranges ob - 
served for nearby SO/a-Sb galaxies by Mollenhoff & Heidtl OOOll) . 
Finally, a background level and photon noise were added to the 
simulated images in order to mimic the instrumental setup and 
signal-to-noise of the WFCAM image. The relative errors in the 
fitted parameters of the artificial galaxies were estimated by com- 
paring the input and output values and were assumed to be nor- 
mally distributed. The standard deviation was adopted as the la 
error in the relevant parameter for the bulge-disc decomposition. 
The best-fitting values and their 3u errors are /i c = 17.07 ± 0.64 
mag arcsec -2 , r c — 11.5 ± 6.0 arcsec, n = 3.46 ± 0.42, 
PA b = 72? 8 ± 0?9, and q b = 0.88 ± 0.03 for the bulge, and 
= 17.06 ± 0.72 mag arcsec" 2 , h = 35.4 ± 18.0 arcsec, 
PA d = 165? 7 ± 6?6, and q d = 0.47 ± 0.06 for the main disc. 

Systematic errors given by a wrong estimation of the PSF 
FWHM and sky level are the most significant contributors to the 
error budget, since the spiral arms do not affect the result and no 
further component is observed. We estimated a 3<r error of 2 and 1 
per cent for the sky level and PSF FWHM, respectively. We anal- 
ysed the artificial galaxies by adopting the correct sky level and a 
PSF FWHM that was 2 per cent larger (or smaller) than the actual 
one or the correct PSF FWHM and a sky level that was 1 per cent 
larger (or smaller) than the actual one. To derive the intrinsic shape 
of the bulge we are interested in q b , qd, PAb, and PAd. Their 3<j 
errors including systematics are smaller than 15 per cent. 



2.2 Intrinsic shape of the bulge 

The method developed by MA+10 was applied to derive the intrin- 
sic shape of the bulge of NGC 4698. It is based upon the geometri- 
cal relationships between the projected and intrinsic shapes of the 
bulge and its surrounding disc. The bulge is assumed to be a triaxial 
ellipsoid with semi-axes of length A and B in the equatorial plane 
and C along the polar axis. The bulge shares the same centre and 
polar axis as its main disc, which is circular and lies on the equato- 
rial plane of the bulge. The intrinsic shape of the bulge is recovered 
from the bulge ellipticity parameter e = (1 — ?b)/(l + ?b) — 0.13, 
the twist angle 5 = 180° - |PA b - PA d | = 87? 1 between the ma- 
jor axes of the bulge and main disc (i.e., the line of nodes), and the 
main disc inclination 9 = arccos qd = 61? 7. 

The relation between the intrinsic and projected parameters 
depends only on the spatial position of the bulge, i.e., on the angle 
^ 4> ^ 90° measured on the bulge equatorial plane between the 
principal axis corresponding to A and the line of nodes. The equa- 
torial ellipticity Z = B 2 /A 2 and intrinsic flattening F = C 2 /^4 2 
of the bulge are given by: 




cos 2 



-F = - sin 4> B cos (20c - 4>b) (1 + Zy 
+ sin (2<f>c - <t> B )sJ(l - Zf - sin 2 (f> B (1 + Z) 2 , (1) 



where: 



B/A 



Figure 2. The distribution of the axial ratios of the bulge of NGC 4698. The 
median value oi B/A and C / A is marked with a segment in the bottom and 
right-hand panel, respectively. 



i e sin 25 ,„„ 

<j) B = arctan —— — = 1.8, (2) 

cos a (1 + e cos 2o) 

1 2 e sin 28 cos 8 „„„„ 

4>c = ■= arctan =- = 89? 6. (3) 

2 e cos 2(5 (1 + cos 2 6) - sin 2 6 

Since the axial ratios B /A and C /A are both functions of the same 
variable 0, their probabilities are identical, i.e., for a given value of 
B/A with probability P(B/A), the corresponding value of C/A 
obtained by Eq. [I]has a probability P(C/A) = P(B/A). This 
allows to obtain the range of possible values of B/A and C/A 
for the bulge of NGC 4698 and to constrain its most probable in- 
trinsic shape by adopting the probabilities P(Z) and P(F) de- 
rived by following MA+10. Fig.f2]shows the distribution of B/A 
and C/A calculated via Monte Carlo simulations. We randomly 
chose 5000 geometric configurations assuming a Gaussian distri- 
bution of gb , qd, PAb, and PAd taking into account their 3<r errors. 
For each geometric configuration we derived 1000 values of B/A 
and C / A according to their probability distribution functions (see 
MA+10). Therefore, the resulting median values of the axial ratios 
and the 3cr confidence intervals from their cumulative distribution 
are B/A = 0.95± ;gf and C/A = 1.60± ;^, respectively. The 
bulge of NGC 4698 is elongated orthogonal to the main disc and 
the probability that it has an equatorial section {B/A < 0.95) is 50 
per cent. 



3 THE NUCLEAR STELLAR DISC 
3.1 Hubble Space Telescope imaging 

The Wide Field Planetary Camera 2 (WFPC2) images of 
NGC 4698 obtained with the filters F450W (Prop. Id. 9042, P.I. 
S. J. Smartt), F606W (Prop. Id. 6359, P.I. M. Stiavelli), and F814W 
(Prop. Id. 9042, P.I. S. J. Smartt) were retrieved from the HST Sci- 
ence Data Archive. The ^606^ images were taken by centring the 
galaxy nucleus on the Planetary Camera (PC), whereas the others 
were taken with the Wide Field Camera (WFC). The PC and WFC 
detectors are Loral CCDs with 800 x 800 pixels and a pixel size of 
15 x 15 ^tm 2 . The PC image scale of 0.046 arcsec pixel -1 yields 
a field of view of about 36 x 36 arcsec 2 . Each WFC detector cov- 
ers 72 x 72 arcsec 2 with 0.091 arcsec pixel -1 . To help in iden- 
tifying and correcting cosmic ray events, different exposures were 
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Table 1. Photometric parameters of the nuclear stellar disc. 
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Figure 3. Isophotal parameters of the nuclear region of NGC 4698 as a 
function of the isophotal semi-major axis based on the analysis of the 
surface-brightness distribution measured in the F450W (left-hand panels), 
F606W (central panels), and F814W (right-hand panels) images, respec- 
tively. From top to bottom: Surface-brightness radial profiles of the galaxy 
(open black squares) and NSD (after convolution with the HST PSF, solid 
blue line). Radial profile of the NSD-to-total luminosity ratio. Radial pro- 
files of the galaxy ellipticity, position angle, and fourth cosine Fourier coef- 
ficient before (open black squares) and after (filled red circles) the subtrac- 
tion of the best-fitting model for the NSD. 



taken with each filter. The total exposure times was 460 s for the 
F450W and F814W filters and 600 s for the F606W filter. The tele- 
scope was always guided in fine lock, giving a typical rms tracking 
error per exposure of 0.005 arcsec. The images were reduced us- 
ing th e CalWFPC reduction pipeline in IRAF (McMaster et al. 
2008). Subsequent analysis including alignment and combination, 
rejection of cosmic rays, and sky subtraction was performed using 
IRAF standard tasks. The flux c alibration to the Vega magnitude 
system was performed following iHoltzman et al. I (fl995h . The flux 
calibration to the Veg a magnitude sys t em in the HST passbands was 
performed following H oltzman et alj i 19951) . 

Following IPizzella et al.l a2002h the unsharp-masked image of 
each WFPC2 frame was built to first gauge the structure and extent 
of the NSD in all the available images. This procedure enhanced 
any surface-brightness fluctuation and non-circular structure ex- 
tending over a spatial region comparable to the standard deviation 
of the smoothing Gaussian. The galaxy nucleus clearly reveals the 
presence of a highly elongated structure in all the images. This is 
the NSD and its location, orientation, and size are the same in all 
the observed passbands. This nuclear structure is associated with 
a central increase in ellipticity, e, and fourth cosine Fourier coef- 
ficient, a4, as measured by performing an isophotal analysis using 
ELLIPSE and shown in Fig. [3] 

The photometric parameters of the NSD were derived in the 
different passbands using the method by_ Scorza & Benderl (1 19951 
hereafter SB95) as implemented by iMorelli et al.l d2004l) although 
with a different treatment of the PSF. The photometric decomposi- 
tion was performed independently for each band-pass. The SB 95 
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method is based on the assumption that the isophotal disciness 
is the result of the superimposition of a host spheroidal compo- 
nent and an inclined exponential disc. Both the components are 
assumed to have perfectly elliptical isophotes with constant but dif- 
ferent ellipticity. The method consists of the iterative subtraction of 
an infinitesimally thin disc model characterised by an exponential 
surface-brightness profile with central surface brightness To.nsd, 
scalelength /insd, axial ratio qnsd, and position angle PAnsd - The 
NSD parameters are adjusted until the departures from perfect el- 
lipses are minimised (i.e., 0,4 ~ 0). To properly derive the pho- 
tometric parameters of the NSD, it is important to account for the 
HST PSF. For each nuclear disc model, the disc-free image of the 
galaxy was obtained from the galaxy image by subtracting the nu- 
clear disc model after convolving with the HST PSF. The adopted 
PSF model was calculated with the TINYTIM package taking into 
accoun t the instrumental setu p and position of the NSD on the given 
image jKrist & Hooklll999l) . The best-fitting value s of the NSD 
param eters and their 3a errors were derived as in IMorelli et al.l 
d2004l) and are listed in Tab.Q] The NSD inclination is calculated as 
insd = arccos qnsd- The NSD is elongated like the galaxy bulge 
as is evident from its position angle. The comparison between the 
isophotal parameters of NGC 4698 measured before and after the 
subt raction of the best-fitt ing model of the NSD are shown in Fig.[3] 
IPizzella et alj fcOO/j) had already analysed the F606W image 
with the SB95 method. They deconvolved the galaxy image from 
the effects of the PSF using the Richardson-Lucy method. We de- 
cided to repeat the photometric decomposition since both an homo- 
geneous photometric decomposition and an estimate of the errors 
on the fitted parameters are required to properly compare the struc- 
tural parameters of the NSD in the different passbands. Moreover, 
this allowed us to verify that the NSD parameters obtained by sub- 
tracting the PSF-convolved disc model from the galaxy image are 
in agreement within the errors with those obtained by subtracting 
the unconvolved disc model from a deconvolved galaxy image. To 
this aim the conversion to Johnson V band has to be taken into ac- 
count for a proper comparison of 7o,nsd and was calculated with 
IRAF task SYNPHOT (V - F606W = 0.27 for a spiral galaxy). 



3.2 Formation of the nuclear stellar disc 

The structural parameters of the NSD (i.e., scalelength, inclination, 
and position angle) are constant within the errors in all the available 
WFPC2 images. Therefore, the location, orientation, and size of the 
NSD do not depend on the observed passband. This also implies the 
absence of colour gradients in the NSD, an important constraint on 
the star formation process. The mean values of the NSD parameters 
are (hnsn) = 48.7 pc, (i NSD > = 76? 1, and (PA NSD > = 69?7. 
The size and luminosit y are consistent with those of the other NSDs 
detected so far (see iLedo et al1l2010L for a census). The NSD is 
oriented as is the bulge, and its major axis is perpendicular to that 
of the main disc, i.e. it is a polar NSD. 

A linear combination of single-age stellar population synthe- 
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sis models was used bv lSarzietaLld2005l) to interpret an HST spec- 
trum measured within the central 0. 13 arcsec (11 pc) of NGC 4698. 
The continuum spectral energy distribution between 3000-5700 
A was best fitted by a very old stellar population (10 Gyr) with 
no evidence for stars younger than 5 Gyr in the galaxy nucleus, 
where the NSD light contamination is large (~ 20 per cent, Fig. 
0. The light fraction from stars younger than 1 Gyr rises to 2 
per cent if supersolar-metallicity models are considered. This im- 
plies that the centre of the NSD is made mostly of old stars. Mod- 
els for chemical and spectro-photometric evolution of galaxy discs 
predict a strong evolution with age of their c olour profiles as di- 
rect consequence of an inside-o ut formation dBoissier & Prantzosl 
1 19991 : IPrantzos & Boissierl boOO ) . In this scenario colour gradients 
are produced early on in the disc centre and propagate outwards. 
The maximum colour gradient measured for the NSD (\d(B — 
R)/dr\ < 0.05 mag kpc -1 ) falls short of the model predictions 
and it is not consistent with a picture wherein the nuclear disc as- 
sembled inside out in the last 5-10 Gyr. 



4 DISCUSSION AND CONCLUSIONS 

The structure of the nucleus and bulge of the Virgo spiral galaxy 
NGC 4698 was investigated through a detailed analysis of HST 
optical and UKIRT near-infrared images, respectively. The galaxy 
is known to host a nuclear disc of gas and stars which is rotat- 
ing perpendicularly with respect to the galaxy main disc. More- 
over, the bulge and m ain disc appear on the sky elongated perpen- 
dicular to each othe r dBertola et aLlll999r. iBertola & Corsinill200d : 
iPizzella et all2002h . 

The equatorial ellipticity and intrinsic flattening of the polar 
bulge were obtained following MA+10 from the apparent elliptic- 
ity of the bulge, twist angle between the bulge and main disc, and 
main disc inclination measured in the UKIDSS near-infrared image 
by adopting the photometric decomposition method by MA+08. 
The bulge of NGC 4698 is remarkably elongated in a perpendicu- 
lar direction with respect to the disc plane (C/A — l.GOtg'}*)- Al- 
though the consistency with axisymmetry was recognised (B/A = 
0.95lQ gj), still the probability that it is significantly triaxial (i.e., 
with an equatorial elliptical section with B/A < 0.95) is 50 per 
cent. The central surface brightness, scalelength, inclination, and 
position angle of the polar NSD in all the available HST images 
were measured by assuming it is an infinitesimally thin exponential 
disc and applying t he photometric deco mposition method of SB95 
as implemented bv lMorelli et al.l d2004h . The size, orientation, and 
location of the polar NSD do not depend on the observed passband, 
as already observed for the few other NSDs fo r which a detailed 
multi-band photometric a nalysis was performed (Kra inovic & Jaffe 
l2004llMorellietal.ll201(j|) . 

The combination of these new results about the complex struc- 
ture of NGC 4698 gives us the opportunity to gain further insight on 
the formation of NSDs with the goal of motivating numerical mod- 
elling to test and refine this scenario. The kinematical decoupling 
between two components of a galaxy suggests the oc currence of 
an accretion event or merging dBertola & Corsinil 19991) . Therefore, 
it is straightforward to explain the existence of the orthogonally- 
rotating dynamically-cold nuclear disc in NGC 4698 as the end re- 
sult of the acquisition of external gas by the pre-existing galaxy. 
Gas dissipation is indeed a necessary ingredient s ince purely stel- 
lar dy namical mergers can not form a nuclear disc (Hartma nn et al.l 
I20T 1). In NGC 4698 the accreted gas settled on the principal plane 
perpendicular to the shortest axis of the triaxial bulge (i.e., perpen- 



dicular to the galaxy main disc) and formed stars. The major axis 
of the NSD is elongated in the same direction as the bulge. The 
link between the angular momentum transport from galactic to nu- 
cle ar scales and formation of NSDs has been recently investigated 
by [Hopkins & O uataertl d2010f). The stellar popu lation in the cen- 
tre of NGC 4698 is very old dSarzi et al.ll2005l) . The absence of 
colour gradients in the NSD is explained if either the star forma- 
tion homogeneously occurred all over the extension of the disc or 
the NSD assembled through an inside-out process that ended more 
than 5 Gyr ago. 
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